Ferromagnetic resonance in systems with competing uniaxial and cubic anisotropics 
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We develop a model for ferromagnetic resonance in systems with competing uniaxial and cubic 
anisotropics. This model applies to i) magnetic materials with both uniaxial and cubic anisotropies, 
and ii) magnetic nanoparticles with effective core and surface anisotropies. We numerically compute 
the resonance frequency as a function of the field and the resonance field as a function of the direction 
of the applied field for an arbitrary ratio of cubic-to-uniaxial anisotropy. We also provide some 
approximate analytical expressions in the case of weak cubic anisotropy. We propose a method 
that uses these expressions for estimating the uniaxial and cubic anisotropy constants, and for 
determining the relative orientation of the cubic anisotropy axes with respect to the crystal principle 
axes. This method is applicable to the analysis of experimental data of resonance type measurements 
for which we give a worked example of an iron thin film with mixed anisotropy. 

PACS numbers: 76.50.+g; 75.75.+a; 75.10.Hk 



I. INTRODUCTION 



The dynamics of magnetic nanoparticles is an area of 
intense theoretical and experimental investigation. From 
the technological point of view, one of the reasons for 
such a great interest partly stems from the growing de- 
mands on the magnetic recording industry. However, 
small nanoparticles, used for high density storage, be- 
come superparamagnetic even at low temperature be- 
cause of the thermal instability of their magnetisation. 
Controlling this effect, in view of room temperature ap- 
plications, requires an understanding of the magnetisa- 
tion dynamics at the nanosecond time scale, taking into 
account the various material properties. The magnetic 
properties of a fine nanoparticle, as compared to the 
bulk material, are dramatically altered due to its reduced 
size and to the boundary effects which are induced by 
the modified atomic environment at its surface. Conse- 
quently, the nanoparticle turns out to be a many-body 
system whose magnetic state cannot, a priori, be faith- 
fully described by a macroscopic approach. On the other 
hand, from the theoretical point of view, investigating the 
dynamic properties of a many-spin particle is a real chal- 
lenge because one is faced with the inherent difficulties 
related to the analysis of the energy potential. This anal- 
ysis is unavoidable since it is a crucial step in the study of 
ferromagnetic resonance (FMR) and the dynamic rever- 
sal of the particle's magnetisation. Nevertheless, there 
exist some regimes of the physical parameters where the 
macro-spin approach may still be used if appropriately 
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extended. For example, in Refs. [l|, 0, Q it has been 
shown, analytically as well as numerically, that when the 
surface anisotropy constant is much smaller than the ex- 
change coupling, the surface anisotropy contribution to 
the particle's energy is of 4*''-order in the net magnetisa- 
tion components and 2"'^-order in the surface anisotropy 
constant. This means that the behaviour of a many- 
spin particle with uniaxial anisotropy in the core and 
relatively weak surface anisotropy (transverse or Neel) 
can be modelled by that of a macro-spin system whose 
effective energy contains an additional cubic-anisotropy 
term. On the other hand, magnetocrystalline anisotropy 
of 4"^-order naturally arises in magnetic materials (e.g., 
iron and YIG) and adds to the 2"'^-order contribution, 
though with an order of magnitude smaller [see Ref. [J]]. 
In iron magnetic multilayers and thin films we also ob- 
serve a mixing of uniaxial and cubic anisotropies which 
arise from their growth on particular substrates, notably 
GaAs (001) [|. 

In Ref. [i*! the effective macro-spin approach was 
used to interpret the /i-SQUID measurements of the ZD 
switching field (or Stoner-Wohlfarth astroid) for a cobalt 
particle of 3 nm diameter (containing 1500 atoms). The 
authors concluded that the cluster-matrix interface may 
be responsible for the main contribution to the magnetic 
anisotropy of the nanocluster. This gives a further in- 
dication that the macro-spin approach with an effective 
potential energy provides a more reasonable approxima- 
tion to the initial many-spin particle than the macro-spin 
Stoner-Wohlfarth model. This may thus be used to in- 
vestigate, though in a phenomenological manner, the dy- 
namics of the particle, and in particular the FMR char- 
acteristics. 

Accordingly, we first consider the general situation 
of either i) a single magnetic moment, regarded as a 
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macro-spin representation of a nanoparticle, with an 
effective potential containing both a 2"^^- and a 4**^- 
order anisotropy terms, including the applied magnetic 
field, or ii) a magnetic material with competing mag- 
netic anisotropics (uniaxial and cubic). We compute the 
resonance field and frequency as functions of the inten- 
sity of the cubic-anisotropy contribution and the static 
magnetic field (with varying direction and magnitude). 
In addition, we also consider the possibility of having 
the cubic-anisotropy axes rotated to some angle with re- 
spect to the axes of the crystal lattice, while the uniaxial 
anisotropy easy axis is maintained parallel to the z axis 
throughout. 

This work is organised as follows: we first outline the 
basic framework of the free energy and parametrise it into 
dimensionless quantities which is convenient for the nu- 
merical calculations. The free energy, containing a uniax- 
ial and cubic anisotropy terms together with the Zeeman 
contribution, is then used to evaluate the ferromagnetic 
resonance conditions. We perform two types of calcula- 
tion: i) the resonance field as a function of the direction 
of the applied field, and ii) the full frequency spectrum as 
a function of the applied field. In each case we have var- 
ied the relative strengths of the anisotropy constants to 
demonstrate how the resonance condition changes with 
the growing influence of the cubic-anisotropy contribu- 
tion. After discussing the results for the general case, we 
consider the specific case of a cobalt nanoparticle as a 
model system, with the parameters of Ref. to esti- 
mate the range of the resonance frequency and field. We 
also estimate the anisotropy constants in an iron thin 
film with both uniaxial and cubic anisotropics by fitting 
experimental data. 



The magnetic moment m experiences the effective field 
defined by Heff = —jj-5E/6s and which is normalized 
with respect to the second-order anisotropy field 



2K2V 



(2) 



leading to the dimensionless field vector 



hoff = /ie,, + /C2(s -11)11 - C sl,R^°'ea, (3) 



where 



i?2 K2' 



(4) 



and the parameter k2 — 0,1 is inserted to allow us to 
switch on or off the uniaxial anisotropy. Using these di- 
mensionless quantities, the energy ([T]) becomes 



Using the parametrisation siO^ip) and e/i(0;i, (p;i)the 
(dimensionless) FMR equation reads [see e.g., Ref. 0] 
and references therein] 



sm 



{dl£) {dl£) - {d^sf 



(6) 



with the right-hand side being evaluated at the equilib- 
rium state So{9o,(po). We have also introduced the re- 
duced angular frequency 



II. ENERGY AND PHYSICAL PARAMETERS 



A. Basic expressions and notation 



with 
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We define the system's net magnetic moment as m = 
/is s with Us = myV, where rriy is the magnetic moment 
density, s (|s| = 1) is the unit vector in the direction of 
m. The energy of m reads 

E = -M.i? (s • e,) - K2V (s -11)2 + :^ (4, + 4, +4), 

(1) 

where n is the uniaxial anisotropy easy axis {K2 > 0) 
and Bfi the unit vector along the applied field. In 
Eq. ([T]) we have used a different form for the cubic 
anisotropy from that often found in the literature, i.e., 
-K4V {slsl + s^sl + slsl), with K4 > 0. These forms 

are related through the identity 1 = {s^ + Sy + s^)^ = 
{s'^Sy + Sysl + s^s^), wherc the first equal- 
ity is due to |sl = 1. The {x'y'z') coordinate system is 
deduced from (xyz) by a rotation of ip around the z axis, 
i.e., Sa' — ^i3=x y z^"^ ^l3^ whcrc R"^ are the matrix 
elements of the corresponding rotation. 



fJ-s 



1 



(8) 



2-fK2V -fH2 
being the scaling time. So, the reduced frequency reads 



UJs_ 

2tt 



2tt 



As mentioned in the introduction, the appearance of 
the 4*''-order term in Eq. ^ may be of two origins: i) 
the natural magnetocrystalline cubic-anisotropy contri- 
bution that adds to the uniaxial anisotropy in real ma- 
terials 0, H, 01 ii) been shown in Refs. [l], [3, H, 
the contribution induced by the surface anisotropy in a 
nanoparticle cut from a cubic lattice. 

While the first situation is quite common in magnetic 
materials, and in particular in thin films, the second 
case deserves further explanation and a review of the re- 
cent results on which the underlying assumption is based. 
This is done in the following section 
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B. Effective energy of a nanoparticle 

Investigating the dynamics of a nanoparticle taking ac- 
count of its intrinsic properties, such as its size and shape, 
crystal structure and surface anisotropy, would require 
the use of an "atomistic" approach with the atomic mag- 
netic moment as its building block. However, within this 
approach one is faced with complex many-body aspects 
with the inherent difficulties related with analysing the 
energyscape (location of the minima, maxima, and sad- 
dle points of the energy potential) . This analysis is un- 
avoidable since it is a crucial step in the calculation of 
the relaxation time and thereby in the study of the mag- 
netization stability against thermally-activated reversal. 
In view of such difficulties, one may then ask the ques- 
tion as to whether there exists an intermediate approach 
with the relative simplicity of the macroscopic approach 
and richness of the many-spin approach, namely a macro- 
scopic model which captures some of the intrinsic features 
of the magnetic nanoparticle. In Refs. [H, 0, S analytical 
as well as numerical calculations showed that a many- 
spin particle, cut from a cubic lattice, when its surface 
anisotropy is small with respect to the exchange cou- 
pling, i.e., when its magnetic state is not far from the 
collinear state, may indeed be modeled by an effective 
one-spin problem (EOSP), i.e., a single macroscopic mag- 
netic moment m representing the net magnetic moment 
of the many-spin particle. The energy of this EOSP (nor- 
malized to JA/", where J is the substance bulk exchange 
coupling and TV the total number of spins in the cluster) 
may be written as 

fEOSP = fc + fl°^+4'^+ff . (9) 

where i) £c is the pure core anisotropy contribution that 
may be uniaxial, cubic, bi-axial, etc. ii) is the pure 
surface contribution that stems from surface anisotropy; 
it is quadratic in the single-site surface anisotropy con- 
stant Kg and quartic in the components of m; it is also 
proportional to a surface integral that depends on the 
size, shape, and crystal structure of the initial many- 
spin particle [see Eq. PTI| below], iii) the contribution 
E^f^ is induced by elongation (or shape anisotropy); this 
term is quadratic in the components of m and linear in 
Ks- iv) arises from a competition between the sur- 
face anisotropy which induces spin disorder that tends to 
propagate deep into the particle, and the core anisotropy 
that tends to expel such spin-noncoUinearities out to the 
particle's border. This core-surface mixing contribution 
is linear in the core anisotropy constant Kc, quadaratic in 
Ks^ and its dependence on m is given by a function that 
mixes quadratic and quartic anisotropics. However, the 
contribution s'^'^ is only relevant if [Kd J) ^ {Ks/ JY- 

Collecting all contributions, a many-spin particle satis- 
fying the above-mentioned conditions, may be described 
by a macroscopic magnetic moment m, representing the 
net magnetic moment of the particle, whose (effective) 
energy may be written as in Eq. with k2 , C being 



regarded as the effective uniaxial and cubic anisotropy 
constants, respectively. It is then important to note that 
the magnitude and sign of these constants depend on the 
intrinsic properties of the initial many-spin particle, such 
as the crystal structural, size and shape, and physical pa- 
rameters such as the single-site surface anisotropy (in in- 
tensity and model). However, we should emphasize that 
these results hold for cubic crystal lattices and quadratic 
surface anisotropy models, such as Neel's or transverse. 

For further reference, we note that in Ref. [H an analyt- 
ical expression was given for the effective constant K^s 
of the surface-induced cubic-anisotropy term s'^^ , when 
the core anisotropy is absent, that is 

K,G = K^, (10) 

where z is the coordination number and n a surface inte- 
gral that depends on the underlying lattice, shape, and 
size of the particle and also on the surface-anisotropy 
model. For a spherical particle (of ^ 1500 spins) cut from 
a simple cubic lattice and with Neel's surface anisotropy, 
K ~ 0.53465. 

Finally, we note that the shape anisotropy is included 
in the uniaxial anisotropy contribution. We also assume 
that in ellipsoidal particles the magneto-crystalline easy 
axis is parallel to the direction of the major axis. 

III. RESULTS AND DISCUSSION 

Before we discuss our results, some remarks are in or- 
der concerning the general numerical method used here. 
In order to obtain, for instance, the resonance field from 
Eq. ([S]), for a given angular frequency ojs, one has to 
find the equilibrium state of the system for a given set 
of the physical parameters k2TC,,h^9h,^h- However, as 
it is not possible to obtain in general analytical expres- 
sions for the equilibrium state, or the global minimum of 
the energy ([5]), we resort to numerical approaches. Ac- 
cordingly, since we only need the absolute minimum, an 
adequate numerical method is provided by the standard 
Metropolis algorithm with random increments, which is 
a global method. Next, once the global minimum has 
been found, we proceed with a fine search by solving the 
Landau-Lifshitz equation with weak damping using the 
Runge-Kutta method. 

In Fig. [l]we plot the resonance frequency as a function 
of the applied field for various values of C, with the cubic- 
anisotropy axes making an angle ip = 0, 7r/4 with respect 
to the (x, y, z) frame. We now discuss some features com- 
mon to both cases oi ip = Q and 7r/4: i) In the case of 
zero field, using the FMR equation ((51), or the effective 
field ([3]), with small C, the global minimum is predomi- 
nantly determined by the uniaxial anisotropy, such that 
cjg ~ fc2 — C- This explains the decrease of Us with C, 
at zero field, since the uniaxial- and cubic- anisotropy 
terms contribute with opposite signs [see Eq. Ob- 
viously, this behaviour does not change when the cubic- 
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FIG. 1: Resonance frequency as a function of the applied field for various values of ( with the cubic-anisotropy axes making 
an angle — 0, 7r/4 with respect to the {x,y,z) frame. 



anisotropy axes are rotated by some angle around the z 
axis, as can be seen in Fig. [1] (right), due to the rota- 
tion symmetry in the equator {xy plane) when h — Q. 
ii) In strong fields, 6q ^ 9h, fo and we can expand 

the FMR equation ([6]) in {(/h), assuming weak cubic 
anisotropy, to obtain the asymptotic behaviour (dashed 
lines in Fig. [1] (left)). This leads to 

h /a \ C , , 

r = I a = cos^ ip + sin^ ^p , b — cos Aijj. 

The coefficient of the term in is {r+l/r)/2 for t/; ~ 
and (r — l/2r)/2 for ip — tt/A. Hence, for -0 = 7r/4 the 
effect of the correction term in C^/h is negligible. Indeed, 
as can be seen in Fig. [T] (right), the high-field asymptote 
is almost independent of C and is approximately given 
by the asymptote for the uniaxial anisotropy, i.e., by 
Eq. (HID with C = 0. 

A point on the curves ujs{h) that is of special impor- 
tance in practice, is that defined by {h = he, tOg — 0) [see 
discussion below] . In Fig. [1] the critical field at which ujs 
vanishes is given by 

/ic = fc2 + C [cos"* -I- sin'' V] • (12) 

In particular, we have he — k2 + C in Fig. [T] (left) and 
he = k2 + C/2 in Fig. [1] (right). This, together with the 
analysis for h — and the asymptote pT|) for strong 
fields, shows that the effect of the ip rotation is to reduce 
the influence of cubic anisotropy. In the case of pure 
cubic anisotropy (^2 = 0) the curves of Ug versus h cross 
the other curves with fc2 7^ 0, and for the sake of clarity 
were not included in Fig. [TJ One should note that in 
fact the critical field h^ at which tOs tends to zero, is the 
field on the Stoner-Wohlfarth curve (or astroid) at which 
the metastable minimum merges with the saddle point 
and loses its local stability. Indeed, from Eq. ^ we see 
that the condition = is just the definition of this 
inflection point. 



The above results can be used to interpret FMR mea- 
surements with sweeping frequency as obtained, for ex- 
ample Jjy the Network Analyser FMR (NA-FMR) tech- 
nique @, Brillouin Light Scattering (BLS) [9] and Pump- 
Probe Microscopy (PPM) [lO|. More precisely, one can 
extract the anisotropy constants K2,K4, and also the 
rotation angle ip. For this, we can use the three inde- 
pendent conditions provided by i) tUg at zero field which 
is ^2 — Ci ii) the high-field asymptote (|lip . and iii) the 
critical point he given by at which a;^ tends to zero. 
Therefore, for a given material with given uniaxial and 
cubic anisotropics, in which the cubic-anisotropy axes are 
at some arbitrary azimuthal angle ip with respect to the 
crystal axes, we can uniquely determine the three pa- 
rameters {K2, K^^^tp). It should be noted, however, that 
these measurements must be made with the field applied 
along one of the hard axes in order to obtain the three 
above conditions. For any other orientation of the ap- 
plied fleld we lose condition iii), because the cusp is no 
longer well defined. In the general situation, where the 
uniaxial-anisotropy axis is tilted at some angle with re- 
spect to the z axis ([001] axis), the above three conditions 
should be rederivcd. 

Fig. [5] shows plots of the resonance field versus the di- 
rection of the applied field for various values of C, and 
for the cubic-anisotropy axes parallel and at ip — 7r/4 
with respect to the {x,y,z) frame. These plots are for 
the fixed frequency LOg = 1.5, which corresponds to 13.5 
GHz in physical units. In effect, this corresponds to cut- 
ting a horizontal line (at w = 1.5) in the frequency- field 
curves (Fig. [T]). The curves with fc2 = 0, ^ = 1 corre- 
spond to the case of pure cubic anisotropy and are in- 
cluded only for comparative purposes. At Oh — 0, we 
have /iros — to — (^2 ~ C)- In strong fields, we can compute 
the resonance field, in principle, by replacing in Eq. (fTT|) 
ujs by a fixed frequency and solving for h. This should 
lead to an approximate expression, albeit somewhat cum- 
bersome, for the resonance field at 6h = tt/2. For this 
orientation, we see that the difference induced by the cu- 
bic anisotropy (C) is suppressed by the ip rotation. As 
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FIG. 2: Resonance field versus the direction of the applied field for various values of ^. The field is applied in the xz plane. The 
frequency is set to uos = 1.5, which corresponds to 13.5 GHz. Also shown for comparison is the case of pure cubic anisotropy 
(^2 = 0,(" = 1). The cubic-anisotropy axes are at i/i = 0,7r/4 with respect to the {x,y,z) frame. (The angle Oh is measured in 
radians.) 



for the frequency plots (Fig. [T]) , comparing the resuhs in 
Fig. [2] we see that the effect of the "0 rotation is again to 
reduce the influence of cubic anisotropy. Angular FMR 
studies in thin films also exhibit the angular variations 
indicated in Fig. [2] [11]. 

By way of illustration and to give some orders of mag- 
nitude of the various constants of the resonance frequency 
and field, we first consider the example of a nanoparti- 
cle, and in particular a (faceted) truncated octahedral 
cobalt particle (of 3 nm in diameter) as obtained experi- 
mentally in Ref. 4. In this reference, the 3D experimen- 
tal switching field was fit to Eq. ^ with an additional 
(relatively weak) anisotropy term along the hard-axis y. 
According to the estimations given in K2 — 2.2 x 10^ 
J/m'^ and K4 ~ 0.1 x 10^ J/m'^, which in our case yields 
C — K4/K2 — 0.05. This implies that the effective 
anisotropy of this nanoparticle is mainly uniaxial. On 
the other hand, using Eq. (jlOp and the parameters given 
thereafter, we estimate the surface anisotropy constant 
as Kg ~ 10~^^ J/atom (or 0.1 erg/cm^), for the above 
mentioned cobalt particle of about 1500 atoms. Similar 
values have been quoted by several authors using Neu- 
tron (quasi)-Inelastic Scattering T5| and FMR on 
cobalt particles. ^From Eq. ((8]), using Hs — 1-4 x 10^ 
A/m we obtain Ts ~ 1.8 x 10"^"'^ s, and hence the angular 
frequency Ug = 1.0 (in Fig. [T]) corresponds to ^ 9 GHz, 
and hj-cs = 1 (in Fig. [2]) to Hicb — /ircs-ff2 — 0.3 T [see 
Eq. ©]. 

At present there are no experimental frequency-field 
data available on nanoparticle systems for comparison 
with theory. However, future work is expected to address 
this issue. 

As a second illustration we show the fitting procedure 
described earlier. For this we have used the frequency- 
field data from Ref. for a 40-monolayer iron film, ob- 
tained by NA-FMR, which exhibit mixed uniaxial and 
cubic anisotropics. This film is part of a bilayer system 
which may have some weak coupling. While the cou- 
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FIG. 3: Comparison of experiment and theory for the res- 
onance frequency as a function of the applied field in a 40- 
monolayer iron thin film, as given in Ref. 



pling itself is important in general, we have neglected its 
effect and only used this data as an illustration of the 
fitting procedure. The results are shown in Fig. [3] where 
full circles represent the experimental data while the line 
corresponds to the high-field expansion of Eq. pT|) . which 
we express as 



ai 02 



(13) 



2--sin2(2V;))C 



ai = 

0.2 



= -- [l + 3sin2(2V')C] , 



1 
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2+-sin2(2^) )C 



after setting k2 — 1- This corresponds to the case with 
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mixed uniaxial and cubic anisotropies. The result of the 
best fit yields H2 ~ -0.175T, C ~ 1.7, and tp ~ 17°. 

It is seen that this approach leads to physically realistic 
estimates of the anisotropy parameters. It is also possi- 
ble to consider the demagnetizing field contribution sep- 
arately from the uniaxial anisotropy, and its value can be 
obtained experimentally as given in Ref. 's'. On the other 
hand, the present approach is applicable only at very 
low (or zero) temperature and as such overestimates the 
anisotropy constants for these measurements which were 
taken from experiments at room temperature. Thermal 
effects on the resonance characteristics can be evaluated 
from the calculation of the transverse ac susceptibility 
as shown in Ref. [3. More generally, one also has to in- 
vestigate the thermally-activated reversal of the magne- 
tization and compute the relaxation rate of a magnetic 
moment with mixed uniaxial and cubic anisotropy and 
in a magnetic field applied at an arbitrary angle with re- 
spect to the uniaxial easy axis. This work is in progress 

a. 

On other hand, in the case of magnetic nanoparticle 
assemblies, one has to take account of the effect of inter- 
particle dipole-dipole interactions (DDI) on the static 
and dynamic properties. In this context, and in the 
macroscopic approach, in Ref. 16, the static behaviour of 
the magnetization, as a function of temperature and ap- 
plied field, was investigated taking account of anisotropy, 
DDI and also of the volume and anisotropy axes dis- 
tributions. In the case of weak DDI, practical approx- 
imate analytical expressions were obtained by perturba- 
tion theory which provides a better approximation than 
the mean-field approach used in Ref. [13 [see also discus- 
sion in Ref. .18j . These expressions involve many lattice 
tensors which account for the effect of the demagnetizing 
field and thus describe the change of magnetization in 
prolate and oblate particle systems. They also show how 
the magnetization deviates from the Langevin law in the 
presence of anisotropy and DDI. Next, in Ref. [l^ the 
effect of DDI on the dynamics of the assembly, and in 
particular on the zero-field-cooled magnetization was in- 
vestigated, and an explanation was given for the change 
of behaviour of the maximum temperature as a function 
of the applied field. It was shown that the transverse 
component of the DDI field creates new saddle points in 
the particle's energy and thereby increases the switch- 
ing rate. In addition, the critical (or activation) volume 
that separates the superparamagnetic from blocked par- 
ticles decreases upon increasing the particle concentra- 
tion. It was found that this volume separates the low- 
field regime dominated by the blocked particles from the 
high-field regime dominated by the superparamagnetic 
ones. As such, as the concentration (or intensity of DDI) 
of the sample is increased, the low- field regime shrinks 
and eventually disappears. 

In connection with the present work, the use of the 
developments in Refs.[la.[l9lis necessary in order to study 
the interplay between the effect of DDI and the intrinsic 
properties of nanoparticles modeled as an EOSP with the 



energy in Eq. (O, and the ensuing effects on the FMR 
characteristics. 



IV. CONCLUSION 

We have studied the ferromagnetic resonance of an 
effective magnetic moment in the general situation of 
an energy potential containing both uniaxial and cubic 
anisotropies. In particular, we have computed the reso- 
nance frequency and field as functions of the applied field 
magnitude and direction. These results can be used in 
interpreting the FMR measurements on magnetic mate- 
rials which exhibit both forms of anisotropy. We have 
provided a simple method for estimating the anisotropy 
constants (uniaxial and cubic) as well as evaluating the 
relative orientation of the cubic-anisotropy axes with re- 
spect to crystalline axes. However, this model assumes 
that the uniaxial anisotropy axis is coincident with one 
of the crystalline axes. The general case is currently be- 
ing investigated and will be published in a future work. 
The frequency-field curves, as illustrated here, can be 
obtained experimentally using techniques such as NA- 
FMR, BLS, and PPM. On the other hand, the present 
FMR analysis may also be used to check whether this 
EOSP approach is a reasonable approximation to real 
magnetic systems and especially to nanoparticles. How- 
ever, it is not obvious how to quantitatively distinguish 
between the cubic anisotropy of magnetocrystalline ori- 
gin and the one that is induced by the surface contribu- 
tion, since these two contributions turn out to be of the 
same order of magnitude. Nevertheless, FMR measure- 
ments on well separated nanoparticles of a very small size 
and grown from a magnetic material of negligible cubic 
contribution could provide us with more precise data for 
this purpose. One possible method would be to study the 
FMR spectra as a function of the particle's size. Indeed, 
changing the size should alter the relative contributions 
of surface and core anisotropies and thus indicate the 
origin of each. 

We have used the model presented here to fit the exper- 
imental data from the high-field branch of the frequency- 
field curve of NA-FMR measurements on an iron thin 
film. We observe good agreement between theory and ex- 
periment which illustrates the application of the model. 

In a recent work we have developed a general the- 
ory for computing the whole spectrum of a many-spin 
particle, including core and surface anisotropy, exchange 
and dipolar interactions, and taking into account the par- 
ticle's shape, size, and the underlying crystal structure. 
We then apply this theory to investigate the effects of 
surface anisotropy on the resonance field and linewidth. 
As has been done in Refs. X, 2^ 2> for the static properties 
(hysteresis cycles, energyscape, and magnetic structure), 
a relationship is being established on the dynamic level, 
and in particular in what concerns the FMR character- 
istics, between the many-spin approach and the effective 
EOSP approach employed here. This should also help 



7 



establish the limit of validity of the EOSP model. 
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